Noise parametric identification and whitening for LIGO 40-meter interferometer data 
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One of the goal of the gravitational data wave analysts 
is the knowledge and the accurate estimation of the noise 
power spectral density of the data taken by the detector, be- 
ing it necessary in the detection algorithms. In this paper we 
showed how it is possible estimate the noise power spectral 
density of gravitational waves detectors using modern para- 
metric techniques and how it is possible whiten the noise data 
before pass it to the algorithms for the gravitational waves 
detection. We report the analysis we made on data taken 
by Caltech 40-meter prototype interferometer to identify its 
noise power spectral density and to whiten the sequence of 
noise. We concentrate our study on data taken in Novem- 
ber 1994, in particular we analyzed two frames of data: the 
18nov94.2.frame and the 19nov94.2.frame. We show that it is 
possible to whiten these data, to a good degree of whiteness, 
using a high order whitening filter. Moreover we can choose to 
whiten only restricted band of frequencies around the region 
we are interested in, obtaining a higher level of whiteness. 
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I. INTRODUCTION 

The building of large interferometers is going to reach 
the phase of data taking: TAMA (Japanese) pj is already 
working; GEO (British/German) will begin to take 
data next year; LIGO (U.S.A.) @[|| in 2002; VIRGO 
(French/Italian) § in 2003 (see e.g. §] @ @ for gen- 
eral and exhaustive description of interferometric Gravi- 
tational Waves detectors). 

The large amount of data produced by gravitational 
wave detectors will be essentially noise and, hopefully, 
buried in noise there will be the signal we are looking 
for. Ground-based interferometric detectors are sensible 
to a broad band frequencies (2-3Hz to more than 1kHz) 
in revealing relative displacement of test masses at the 
near and far extremities of interferometer arms due to 
GW signal, but unfortunately a lot of other factors can 
cause a similar displacement. The test masses are sus- 
pended to pendular structures in order to isolate them 



from seismic noise Q9J, but thermal noise of the suspen- 
sion chain will cause a displacement of the mass [ ^0[ . 
Also shot noise and radiation pressure of the laser will 
move the mirrors §] (ll). The physicists are working in 
modeling all possible causes of noise in the interferometer 
giving out a sensitivity curve of the apparatus ]l2| ] |l~3j ] 
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This curve is limited at low frequencies by seismic 
noise; in the middle band by thermal noise and at high 
frequencies (higher than 0.7 — 1kHz) by shot-noise. The 
sensitivity curve of these detector is a broad-band noise 
plus several very narrow peaks due to the violin modes of 
the suspensions wires, that will make the detection of a 
gravitational signal in this frequency band very difficult. 
For this reason efforts have been made in the preparation 
of the analysis of data for cutting |l9| (2(]] out these 
resonances. 

The analysis of data to detect the gravitational sig- 
nal requires an accurate knowledge of the noise, which 
means a statistical characterization of the stochastic pro- 
cess, and in the case of local stationarity and Gaussian 
nature an accurate estimation of the Power Spectral Den- 
sity (PSD). 

Moreover if the output of the interferometer will be non 
stationary over a long period of time, we must be able in 
following the changes in the PSD. A way to achieve this 
is to estimate the PSD on a chunk of data at different 
interval of time, using classical techniques [^4| [p5j p(] ]. 

We proposed the use of adaptive methods to follow on 
line the change in the feature of the spectrum in such a 
way to have at any desired instant the correct curve for 
the PSD g). 

If we are able in identifying the noise of our detector 
we can also apply the whitening procedure of the data. 

The goal of a whitening procedure is to make the se- 
quence of data delta-correlated, removing all the correla- 
tion of the noise. Most of the theory of detection is in the 
frame of a wide-sense stationary Gaussian white noise, 
but in our problem the noise is surely a colored one and, 
in principle, there could be present non stationary and 
non Gaussian features. If we whiten the data, supposing 



1 



hence to be in the frame of a stationary and Gaussian 
noise, we can apply the optimal algorithm detection |p8| . 

Moreover when we are searching a transient signal of 
unknown form it is very important to have a whitened 
noise p9[ . The importance of whitening data is also 
linked to the possibility of reducing their dynamical 
range fj) @. 

In this paper we show how to identify and how to 
whiten the noise data produced by LIGO-40 meter in- 
terferometric detector before applying any algorithm de- 
tection. 

The data taken from Caltech 40-meter prototype in- 
terferometer in November 1994 J2^] have been given 
from the LIGO collaboration at disposal of the data anal- 
ysis groups of other experiments to perform some test of 
data analysis algorithms. A typical PSD of the output 
of LIGO-40 meter data is plotted in figure 111. It is char- 
acterized by a huge number of spectral lines. Part of 
them are due to the violin modes, part to the laser noise 
and mostly of them to the 60Hz and its harmonics power 
supply. 

We use these data to check the algorithms of paramet- 
ric noise identification and the algorithms which perform 
the whitening on realistic data taken by an interferomet- 
ric gravitational waves detector, having already analyzed 
and tuned these techniques on Virgo-like simulated data 
(see reference ). 

The concepts and the notation in this paper are de- 
scribed in detail in our earlier work p7{ . Sections II and 
III give a brief summary of the methods, but the reader 
is directed to Ref. |27j for details. 
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FIG. 1. Example of PSD of LIGO 40-meter data 



II. PARAMETRIC MODELING AND 
WHITENING 

The advantages of power spectral parametric estima- 
tion with respect to the classical spectral methods are de- 
scribed in an exhaustive way in reference p4| . We want 
to underline that one of its advantage in linked to the 
better spectral resolution we obtain if we suppose that 
the process we are analyzing is governed by a dynamical 
law, because we can use the estimation of autocorrela- 
tion function until a certain lag and then extrapolate its 
value at successive lags, under the dynamical hypothesis 
we made. 

An other advantage is that we may compress the infor- 
mation of the PSD in a restricted number of parameters 
(the parameters of our model) and not in the full auto- 
correlation function. This can help us, for example, if we 
want to create a data base of noise sources. 

In this paper we analyze only the Auto-Regressive 
(AR) parametric modeling because, as a byproduct, it 
offers the possibility to write down a linear whitening 
filter |3 ■ 

A. AR models 

An Auto- Regressive process x[n] of order P (AR(P) ) 
is governed by the relation 

p 

x[n] = — a[fc]x[n — k] + w[n] , (1) 

k=l 

which links the value of the data x[n] to the previous P 
ones by the AR parameters a[k], being w[n] a Gaussian 
white noise process. The theoretical form of the PSD for 
such a process is given by 

R ^ |1 + J2k=i afccxp(-z27rfc/)| 2 ' 

being a one of the parameters to be estimated linked to 
the RMS of the process. 

The relationship between the parameters a[k] of the 
AR model and the autocorrelation function r xx (n) of the 
process x[n] is given by the Yule- Walker equations 

\U — f ~ SLi a i r xx[k - I] for k > 1 ,„, 
r xx [fcj - | _ e p =i ^ H] +(j2 forfc = 0i W 

The problem of determining the AR parameters is 
the same of that of finding the optimal "weights vec- 
tor" w = Wk, for k = 1, ...P for the problem of linear 
prediction p4j . In the linear prediction we would pre- 
dict the sample x[n] using the P previous observed data 
x[n] = {x[n — 1], x[n — 2] . . . x[n — P]} building the esti- 
mate as a transversal filter: 
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)[n] = Wf;x[n — k] 



(4) 



k=l 



We choose the coefficients of the linear predictor by 
minimizing a cost function that is the mean squares error 
e = £[e[n} 2 ], being 



en 



x n\ 



(5) 



the error we make in this prediction, obtaining the so 
called Normal or Wiener- Hopf equations 



= r xx [0] - ^2 Wkr X x[-k] , 



k=l 



which are identical to the Yule- Walker ones with 
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(6) 



(7) 

(8) 



This relationship between AR model and linear pre- 
diction assures us to obtain a filter which is stable and 
causal pi| ]: this is the fundamental relation in building 
our whitening filter. 

A method of solving the Yule-Walker equation is the 
Durbin algorithm p^j . 

The strategy of this method is the knowledge of the 
optimal (P — l)th order filter to calculate from it the 
optimal Pth order. 



B. Link between AR model and whitening filter 

When we find the AR P parameters that fit a PSD of 
a noise process, what we are doing is to find the optimal 
vector of weights that let us reproduce the process at the 
time n knowing the process at the P previous time. All 
the methods that involves this estimation try to make the 
error signal (see equation (|^) ) a white process in such 
a way to throw out all the correlation between the data 
(which we use for the estimation of the parameters). So 
one of the output of the filter which identifies the noise 
will be the whitened part of the data. 

The Durbin algorithm, used to estimated the parame- 
ters a[k], introduces in a natural way the Lattice struc- 
ture for the whitening filter (32). The equations which 
represent our lattice filter in the time domain could be 
written 



e f p [n] 



e b p [n]=e b p _ x [n-l] 



1], 



(9) 
(10) 



being e b p the backward error, that is the error we make, in 
a backward way, in the prediction of the data x[n — p+l] 
using p — 1 successive data, and the forward error, 
that is the error we made in estimating the data x[n] 
using the P previous data. The coefficients k p are the 



so-called reflection or PARCOR coefficients, being linked 
to the correlation that the data x[n] has with the x[n— p], 
neglecting the p — 1 values of x in between. In figure || is 
showed how the lattice structure is used to estimate the 



forward and backward errors (see reference \27\) 
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FIG. 2. Lattice structure for Durbin filter. 

Using a lattice structure we can implement the whiten- 
ing filter following these steps: 

• estimate the values of the autocorrelation function 
r ra [fe], < k < P of our process x[n); 

• use the Durbin algorithm to find the reflection co- 
efficients k p , 1 < p < P; 

• implementation of the lattice filter with these coef- 
ficients k p initiating the filter eg[n] = e^n] = x[n]. 

In this way the forward error at the stage P-th is equiv- 
alent to the forward error of a transversal filter and rep- 
resents the output of the whitening filter. 

The procedure of whitening will be accomplished be- 
fore applying the algorithms for the detection of gravita- 
tional signal of different wave forms. The level of white- 
ness of the data needed for the various algorithms could 
be different. We want now to introduce a parameter that 
let us quantify the level of whiteness of data at the output 
of whitening filter. 



C. The "whiteness" of data: measure of flatness of 
PSD 

The spectral flatness measure for a PSD is defined as 



exp(^/^ 2 /2 ln(P(/))rf/) 



(ii) 
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where the integral is extended in the bandwidth of 
Nyquist frequency; this parameter satisfies 



III. ADAPTIVE METHODS:LEAST SQUARES 
LATTICE 



< £ < 1 



(12) 



If P(f) is very peaky, then £ ~ 0, if P(f) is flat than 
£=!• 

Using the definition ( |ll| ) the flatness for a process at 
the output of a whitening filter built with a minimum 
phase filter (as the AR filter is) is 



£e=£- 



[0] 



(13) 



where r xx [0] and r ee [0] are the values of the autocorrela- 
tion function of the process before and after the whiten- 
ing procedure, and £ is the value of flatness for the initial 
sequence pjj. 



D. Order Selection 

The idea of the whitening filter is that the process we 
analyze is an autoregressive one and that once we have 
the AR parameters we can use them in the whitening 
filter. 

In general we don't know the order of our process, even 
if we suppose that it is an AR one. If it is an AR of order 
P, and we use an order p < P, the fitted spectrum will 
be smoother than the original one; if we choose an order 
p > P, there may be spurious peaks in the spectrum. In 
both cases the whitening will be not good. 

If our process is not AR, the number of parameters 
could be in principle infinite. We must then fix a criterion 
that let us select the right order of the process, or at least 
the best one. 

We used the classical order selection criteria p4j pj| 
p3| , that is the Akaike information criterion (AlCj, the 
forward prediction error (FPE) the Parzen's criterion 
(CAT) and the minimum description length (MDL) one 



The implementation of an adaptive filter follows two 
steps: the filtering of the input data and the adjustment 
of the filter parameters with which we process the data to 
the next iteration. The filters parameters are updated by 
minimizing a cost function. The way in which we build 
this cost function distinguishes the adaptive methods |3^] 
|M The Least Squares based methods build their cost 
function using all the information contained in the error 
function at each step, writing it as the sum of the error 
at each step up to the iteration n : 



(18) 



being 



e(i\n) = d[i] 



N 



k=l 



Xi-kWk n 



(19) 



AIC(P) = Aloge(P) + 2P , 



(14) 



where d is the signal to be estimated, x are the data of 
the process and w the weights of the filter. We introduced 
the forgetting factor A that let us tune the learning rate 
of the algorithm. This coefficient can help when there are 
non stationary data in the sequence and we want that the 
algorithm have a short memory. If we have stationary 
data we fix A = 1 . 

There are two ways to implement the Least Squares 
methods for the spectral estimation: in a recursive way 
(Recursive Least Squares or Kalman Filters) or in a Lat- 
tice Filters using fast techniques |27j J3^] . The first kind 
of algorithm, examined in p5[ , has a computational cost 
proportional to the square of the order of filter, while the 
cost of the second one is linear in the order P. We will re- 
port only the application of Least Squares Lattice (LSL) 
algorithm on LIGO 40-meters data, because in our pre- 
vious work |27) the authors showed that this is the best 
choice among the adaptive algorithms the authors ana- 
lyzed for the problem we have to face. 



FPE(P) = e(P) 



N + P+ 1 
N-P-l ' 



(15) 



IV. RESULTS ON 40-METER LIGO DATA 
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MDL(P) = Aloge(P) + Plog N , 



(16) 
(17) 



where e(P) is the mean square error at the order P and 
N is the length of data. In literature the MDL criterion 
is considered the best among them, because it is robust 
with respect to the len gth of the sequence, while the oth- 
ers depend a lot on N [Ea . 



We use two frames of the LIGO 40-meter interferome- 
ter data: the 18nov.2.frame and the 19nov.2.frame. We 
use the GRASP interface to read the data using only the 
frames taken with the interferometer locked j3(]] 36 1 . We 
are not interested to the absolute sensitivity of such inter- 
ferometer, but only in the features of the Power Spectral 
Density, so we don't use calibrated data, and we analyzed 
the IFO output of the data in ADC counts. This is a fast 
channel which was sampled at about 10kHz. 

In figure [l] a typical power spectral density of these 
data is displayed: it is characterized by a high number of 
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lines, which could be mainly harmonics of the electromag- 
netic interference at 60Hz, by two broad peaks around 
180Hz and 300Hz which could be due to the laser power 
supply, by the violin resonances around 500-600Hz, and 
a continuous background, due to electronic noise, shot 
noise, seismic noise etc. 

It is evident that we would need a high number of 
parameter to fit the spectrum and to build the respective 
whitening filter. 

A. 18nov94.2.frame 

We try to obtain an estimate of the order we need 
for whitening the noise PSD of this frame. The MDL 
criterion gives an order of 515. As first test we use this 
value to whiten the 18nov94.2.frame data. 

In figure || we plotted the PSD of 18nov94.2.frame av- 
eraged on 30 sequence of data of ~ 6s in input to the 
Durbin and LSL whitening filters. In the same figure we 
can see the AR(515) fit to the PSD and the PSD at the 
output of the whitening filters. 

In the first figure it is evident that the filters succeed 
in whitening the overall behavior of the noise power spec- 
trum, but not in eliminating all the spectral lines. 

Besides the several harmonics spectral lines, the main 
problem in the whitening of these data is linked to the 
peaks in the range of frequencies from to 1000Hz. These 
peaks are broad, and it is difficult to fit them with an all- 
poles models as it is evident in the zoom we made in this 
band of frequencies in figure 0. 
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FIG. 3. Application of AR(515) whitening filter 



Notwithstanding these problems, the values of flatness 
at the output of whitening filters is quite good as we see 



in figure |^ and in table |. 
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FIG. 4. Flatness at the output of Durbin and LSL whiten- 
ing filter of order P=515. 



18nov94.2 


Durbin 


LSL 


0.0220 


0.877 


0.917 



TABLE I. Values of flatness on averaged psd before and 
after the whitening filters. 
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This is reflected in the time-domain output of the filter 
where it is evident that the rms of the data is remarkably 
reduced (see figure || ) . 
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FIG. 5. Time-domain input and outputs of whitening filters. 

It is to be noted that the results of the adaptive al- 
gorithm reported in figure |J and in table | are a little 
better that those of the 'static' Durbin algorithm. These 
can be explained by the fact that the estimation of pa- 
rameter of Durbin whitening needs a previous estimation 
of autocorrelation function, while the adaptive LSL filter 
estimate the parameters directly by the time data, avoid- 
ing the problem of errors propagation. Moreover in those 
runs the adaptive algorithm was tuned in such a way to 
have a long memory to let the filter continue to adjust 
its parameters with the input of new data. 



1. Whitening filter of higher order 
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FIG. 6. Whitening of 18nov94.2.frame data with a filter of 
order P=662. PSD averaged on 30 sample of 6 sec. 

We try the whitening procedure using firstly a filter of 
P = 662 coefficients and then a filter of P=2682. 

In figure |^ we report the PSD at the output of whiten- 
ing filter of order P — 662 as averaged on 30 samples of 
6s of data. As it is evident there exists always the prob- 
lem in the low frequency part of the spectrum due to the 
two broad peaks. A typical value of the flatness on one 
sample of data is ^0.5 for the Durbin filter and ^0.55 
for the LSL one. 

In figure ^ we plotted the outputs of whitening filters 
with P = 2682 as suggested from the AIC, FPE and CAT 
criteria. The higher order filters succeeded in eliminating 
the narrow spectral lines at higher frequencies, but the 
problems in the band under 1000Hz still remain. The 
flatness of the output PSD is remarkably better in this 
case. A typical value for the Durbin whitening filter is 
£ ~ 0.6, while for the LSL one is £ ~ 0.7 



Until now we used the order of whitening filter as the 
one estimate on the averaged PSD or autocorrelation 
function. Now we run the order selection criteria on a 
single sequence of data long ~ 6 s. The results are re- 
ported in table [n]. 



MDL 


AIC 


FPE 


CAT 


662 


2682 


2682 


2682 



TABLE II. Minimum of order selection criteria on a single 
sequence of data for the 18nov94.2.frame. 
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FIG. 7. Whitening of the 18nov94.2.frame with 2682 pa- 
rameters. 



2. Whitening the data on a restricted band of frequencies. 



We find again that the MDL gives the minimum value 
for P . We whiten the data both with the value found by 
the MDL criterium, that is P = 329, and with the value 
P = 1154 which is the value found by the AIC and FPE 
criteria. 
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By considering the fact that the signal detection algo- 
rithm will act on restricted band of frequencies, where it 
is most probable to find a gravitational signals and that 
the whitening procedure will be limited to an interval of 
frequencies and not on all the Nyquist interval, we try 
to whiten the 18nov94.2.frame on the band of frequen- 
cies from to 1234Hz, where the algorithms have met 
with major difficulties in whitening the noise spectrum. 
To decimate the set of data at our disposal we let the 
data pass through a Butterworth low pass filter of 4 or- 
der that cuts off all the frequencies over 1234Hz and then 
we re-sample the data with a sampling factor 4. 

The output of this pre- filter procedure has been sent to 
the whitening filter. We estimate the best order for the 
whitening filter using a single realization of 6s of data and 
applying to it the order selection criteria (see table [II). 
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FIG. 8. PSD at the output of whitening filters with 
P = 329 for 18nov94.2.frame data in the band - 1233.5Hz. 

The PSD at the outputs of whitening filters with P = 
329 are plotted in figure |[ 

Even if we have restricted the band of interest the dif- 
ficulties of whitening the broad peaks remain. Moreover 
the band of frequencies around 600Hz is characterized by 
several spectral lines plus some broad peaks, due to violin 
resonances jl9| . Most of the parameters of the whitening 
fit try to cancel out the broad peaks, so the whitening 
is worse also on narrow lines. Perhaps it could be useful 
separate the whitening in two steps, eliminate the nar- 
row lines and whiten the broad peaks or augmenting the 
number of parameters. In fact if we use P = 1154 the 
whitening we obtain in this band of frequencies is good 
as you can see in figure ^. 



TABLE III. Order selection criteria for 18nov94.2.frame on 
the band - 1233.5Hz. 
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FIG. 9. PSD at the output of whitening filters with 
P = 1154 for 18nov94.2.frame data in the band - 1233.5Hz. 



FIG. 10. PSD at the output of whitening filters with 
P = 3000 for 18nov94.2.frame data in the band - 1233.5Hz. 



The typical value of the flatness for the whitening with 
P = 1154 is between 0.6 and 0.8 and the Durbin filter 
and LSL filter have similar performance. 

Moreover the choice of order of the filter is not tied up 
by the order selection criteria. 

As we have already outlined p7[ this number is an 
indicative one; the selection criteria depend also on the 
length of sample to analyze. If we use for example frame 
of data long 26s, we find P — 779 for the MDL and 
P = 2694 for the FPE, AIC and CAT. 

Just as a trial we make a whitening with P = 3000 on 
a single set of these data. The result is plotted in figure 
|i"o| . The whitening is very good: in fact the parameter of 
flatness is 0.77 for the Durbin whitening algorithm and 
0.83 for the LSL whitening algorithm. 



This is a proof that we can succeed in whitening the 
data if we choose a sufficiently high number of parameters 
to compensate the fact that we have all-poles models and 
that we need a high number of poles to whiten broad 
peaks. 



B. 19nov94.2.frame 

We perform the same kind of analysis even on the data 
set 19nov94.2.frame. 

We find the order of the filter using the results of the 
order selection criteria on a set of data long ~ 6s, choos- 
ing P = 663 that is the best order find by the MDL 
criterion (see table tV). 



MDL 


AIC 


FPE 


CAT 


663 


2351 


2351 


2351 



TABLE IV. Order selection criteria for 19nov94.2.frame on 
the Nyquist interval. 
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In figure [11] we showed the results of whitening filters 
on 19nov94.2.frame data at an order P = 663. 



In figure |l2| we report the whitening results obtained 
with filter of P = 2351 order. The whitening is remark- 
ably better .The values of the flatness are ~0.6 for the 
Durbin filter and ^0.63 for the LSL one. 



1(TF 



— LIGOPSD 

— Duibin whitening 
— LSL whitening 




1000 2000 3000 4000 

HI I/, i 



FIG. 11. 19nov94.2.frame: psd of input and output for 
whitening filters with P = 663. 

The kind of problem we met in the whitening results 
is the same as in 18nov94.2.frame, that is besides the 
several lines at high frequencies, it is difficult to whiten 
the low frequency broad peaks. A typical value of the 
flatness at the output of Durbin filter is 0.53, whilst for 
the LSL filter is 0.55. 



1. Whitening on a restricted band of frequencies 

We perform the whitening procedure also on this frame 
using a restricted band of frequencies. The selection or- 
der criteria give in this case P = 412 for the MDL and 
P = 912 for the others (see table 0). 



o 10 

K 
I 

o' 




100 200 300 400 500 600 700 800 900 1000 1100 1200 




mi 



FIG. 13. 19nov94.2.frame: PSD of input and output for 
whitening filters on the band — 1233.5 Hz with P = 4f2. 
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CAT 


412 


912 
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TABLE V. Order sefection criteria for f 9nov94.2.frame on 
the band - f 233.5 Hz. 
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FIG. 12. f9nov94.2.frame: PSD of input and output for 
whitening fitters with P = 235 f. 
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FIG. 15. 19nov94.2.frame: psd of input and output for 
whitening filters on the band - 1233.5 Hz with P = 4096. 
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FIG. 14. 19nov94.2.frame: PSD of input and output for 
whitening filters on the band — 1233.5 Hz with P = 912. 

The whitening filter with these order have been im- 
plemented on 30 samples of length 6s. The results are 
plotted in figures ^3||lJ. 

The whitening with P — 912 is good; in this case the 
typical value for the flatness parameter is 0.65 for Durbin 
filter and 0.7 for LSL one. We make the whitening, just 
as a trial, on this set of data using filter of P = 4096. In 
figure [l5] we plotted the results. The whitening procedure 
is good and the values of flatness is 0.66 for Durbin filter 
and 0.77 for LSL one. These values are not very different 
from the ones obtained with filters of order P = 912. 




V. CONCLUSION 

We showed that we can obtain a whitened PSD on re- 
alistic interferometric data. In the two frames 18nov94.2 
and 19nov94.2 of LIGO 40-meter interferometer we pro- 
cessed we found that, restricting our analysis to a band 
of frequencies and using a sufficiently high number of pa- 
rameters, we can obtain whitened data in time domain 
on which we can apply our detection algorithms. It could 
be useful to statistically characterize the LIGO data to 
check out the presence of non stationary or non lineari- 
ties in order to study the behavior of the whitening filters. 
Indeed we verified that the adaptive algorithm has a bet- 
ter performance with respect to the static algorithm so 
it could be useful if we have to face with non stationary 
data. 

In whitening the LIGO 40-meter data we applied the 
filters directly to the data but, in principle, it could be 
useful before applying the whitening filters to cancel out 
the harmonics of 60Hz |o[ @ @- We want to under- 
line that we apply these techniques to check their capabil- 
ity in whitening the power spectral density, but we could 
apply them in a more efficient way if we know in advance 
the kind of analysis we want to perform after the whiten- 
ing, because we can choose to have a smoother level of 
whiteness or to analyze a narrow band of frequencies. 
The filters we use belong to the family of Kalman filter, 
but the dynamical law we suppose to underline the pro- 
cess is a general one, so in principle we can whiten all the 
spectral lines of the PSD. This could be advantageous if 
we want to perform the whitening in one step, because 
we know that the 'risk' of finding a periodic gravitational 
signal with an amplitude comparable to the ones of the 
violin modes or other spectral lines in the PSD is very 
low. 
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